Clinical Pharmacology & Therapeutics

Background: Mocravimod is an oral sphingosine-1-phosphate (S1P) receptor modulator. This Phase 1 multiple-ascending-dose study evaluated its safety, tolerability, pharmacokinetics (PK), and pharmacodynamics (PD) in healthy volunteers. Methods: In this double-blind, randomized, placebo-controlled, parallel-group trial, 60 healthy male volunteers were enrolled in five cohorts. Mocravimod was administered once daily at 0.3, 0.6, 1.2, or 3.0 mg for 14 days, or at 2.0 mg for 28 days. Safety assessments included adverse events (AEs), laboratory tests, vital signs, electrocardiography, and Holter monitoring. PK of mocravimod and its active metabolite, mocravimod-phosphate, and PD effects on absolute lymphocyte count (ALC) and leukocyte subsets were assessed. Results: Fifty-nine of 60 participants completed the study. One participant in the 3.0 mg cohort discontinued treatment because of asymptomatic transaminase elevation. No deaths or serious AEs occurred. AEs were mostly mild or moderate, transient, and showed no clear dose relationship. Mocravimod produced dose-dependent reductions in ALC from 0.6 mg onward, with maximum geometric mean reductions of 65%, 74%, 83%, and 77% at 0.6, 1.2, 2.0, and 3.0 mg, respectively. ALC values recovered to above the lower limit of normal during follow-up in all cohorts. Holter monitoring showed an initial placebo-corrected reduction in heart rate of approximately 10-15 beats/min at doses of 1.2-3.0 mg, which attenuated with continued dosing. One participant in the 3.0 mg cohort had a recurrent daytime second-degree atrioventricular block (Mobitz I/Wenckebach), reported as a mild non-dose-limiting AE. No QT prolongation was observed. Exposure to mocravimod and mocravimod-phosphate increased approximately dose-proportionally. Steady state was reached by Day 14 (Day 28 in the 2.0 mg cohort), accumulation was approximately five- to sevenfold, terminal half-lives were approximately 100-40 hours for both analytes, and parent-to-metabolite exposure ratios were close to 1. Conclusions: Once-daily mocravimod up to 3.0 mg for 14 days and 2.0 mg for 28 days was generally well tolerated and showed predictable S1P-modulator class effects on lymphocyte counts and heart rate, with PK properties supporting once-daily dosing and further clinical development.

Background Classical pharmacokinetic-pharmacodynamic (PK/PD) theory models exposure-effect in two dimensions: magnitude and time. Rate-dependent toxicity has been documented across therapeutic domains but never formalised as a conserved biological constraint. Methods We developed the Human Adaptive Rate Limit (HARL) framework, formalising the maximum tolerable velocity as |dS/dt|_max = sigma_max / tau. We validated HARL across five domains using published trial data and a reanalysis of the longitudinal biomarker data from the 202-patient CAR-T cohort of Wei et al (2023). An 8-ODE quantitative systems pharmacology model guided biomarker selection. Early biomarker velocities (maximum positive slope within days 0-5) were computed for ferritin and D-dimer. Patients were classified as high-risk only if both velocities exceeded their thresholds (dual-velocity classifier). Thresholds were identified by grid-search optimisation of the Youden index and assessed by leave-one-out cross-validation. Findings A prospective crossover study (Kleinbloesem 1987, n=8) demonstrated that matched steady-state nifedipine concentrations produce divergent haemodynamic responses depending solely on rate of rise, anticipating the dose-related mortality signal subsequently reported across ~8350 patients with coronary heart disease (Furberg 1995), a meta-analysis that was itself debated. Convergent evidence spans haematology (CHOIR, 1432 patients, hazard ratio [HR] 1.34 [1.03-1.74] for aggressive Hb correction), radiation (dose-rate effectiveness factor [DDREF] 1.5-2.0), and infusion pharmacology. In the CAR-T cohort, high-risk classification (ferritin >232 ng/mL per day AND D-dimer >1.21 mg/L per day) predicted severe CRS with 100% sensitivity (~78% specificity) in safety rule-out mode and 91.1% sensitivity (93.6% specificity, AUC 0.95 [95% CI 0.91-0.98]) in Youden-optimised mode. Median kinetic lead time was 4 days (range 3-7) before clinical decompensation. Interpretation Biological tolerability is three-dimensional. HARL unifies rate-dependent toxicity across domains spanning minutes to weeks. MTDyn--specifying target level and allowable rate of change--should supplement conventional dose-response assessment.

Background: This study aimed to evaluate real-world adverse event (AE) signals of EV to provide evidence-based guidance for its safe clinical application. Methods: Data from the FDA Adverse Event Reporting System (FAERS) database from the period of 2019 Q1-2025 Q3 were analyzed. Disproportionality analysis algorithms, including the reporting odds ratio (ROR), proportional reporting ratio (PRR), Bayesian confidence propagation neural network (BCPNN), and empirical Bayes geometric mean (EBGM), were utilized to mine safety signals.The time to onset (TTO) was evaluated using the Weibull distribution model. Results: Among 11,697,906 reports, 4,177 EV-treated patients experienced 14,511 AEs. The most common System Organ Classes (SOCs) were skin and subcutaneous tissue disorders (18.23%), general disorders and administration site conditions (13.17%).Multi-algorithm consensus identified 179 positive signals. Alongside known toxicities (rash, peripheral neuropathy, hyperglycemia), potential new signals emerged, including dysgeusia, atypical skin lesions, and myelosuppression. Median TTO was 14 days, with the Weibull {beta} of 0.736, confirming an "early failure" profile. Subgroup analysis revealed toxicity heterogeneity: patients aged [&ge;]65 and females exhibited stronger signals for fatal severe cutaneous adverse reactions, while patients aged < 65 and males showed higher susceptibility to neurological and metabolic toxicities. Conclusions: The real-world safety profile of EV confirms known toxicities, reveals new risks (e.g., dysgeusia), and shows toxicity concentrated in the first treatment cycle. Clinical practice requires proactive monitoring during the first two weeks using demographic-specific strategies: vigilance for fatal skin toxicity in elderly and female patients, and close follow-up of neurological and metabolic indicators in younger and male populations.

Abstract Introduction The sharing of individual patient data is essential for advancing pharmacometrics but is strictly limited by privacy regulations (e.g., GDPR). While synthetic data generation offers a legally compliant alternative, its structural impact on complex nonlinear mixed-effects (NLME) modelling remains largely unexplored. This study aimed to benchmark five generative artificial intelligence algorithms by evaluating the balance between data privacy and the preservation of structural PK properties and clinical dosing guidance. Material & methods A daptomycin two-compartment PopPK model was used to simulate a reference cohort of 500 patients. Five generative algorithms (Modified AVATAR, Gaussian Copula, Synthpop, TVAE, and CTGAN) produced 100 independent synthetic datasets each. A two-stage evaluation framework was applied: first, a statistical indistinguishability test based on logistic regression (AUC ROC) was used as a macroscopic pre-selection criterion to determine algorithm eligibility for NLME modelling and privacy risk assessment. Privacy risk was independently quantified using the Anonymeter framework (Singling Out and Linkability attacks). Eligible algorithms were further evaluated on PK parameter recovery bias and clinical dosing simulations. Results Deep learning architectures (TVAE, CTGAN) were excluded at the pre-selection stage due to both biologically implausible covariate generation and high macroscopic detectability (mean AUC ROC = 0.837 and 0.986, respectively). Synthpop, AVATAR, and Gaussian Copula all passed the indistinguishability threshold (AUC ROC = 0.475 +- 0.033, 0.490 +- 0.013, and 0.619 +- 0.031, respectively) and proceeded to NLME evaluation. However, attack-based privacy assessment revealed that Synthpop carried an unacceptable singling-out risk (0.035), disqualifying it from privacy-preserving data sharing. AVATAR and Gaussian Copula demonstrated acceptable privacy profiles (singling-out = 0.004 and 0.001; linkability = 0.010 and 0.003, respectively). At the structural level, Gaussian Copula injected stochastic noise inflating residual error (+157.0%) and V1; (+25.9%), blunting predicted Cmax and predisposing to empirical dose escalation and risk of toxicity. AVATAR acted aSs a smoothing filter, deflating V2; (-48.3%) and underestimating CL (-12.9%). Forward clinical simulations confirmed directionally opposed prediction errors: Gaussian Copula consistently underestimated Cmax across standard and renally impaired profiles (-14.5% and -16.0%, respectively), predisposing to empirical dose escalation, whereas AVATAR- and Synthpop-derived models overestimated Cmax and Cmin in the obese infected patient (+14.7% and +8.2%, respectively), compounding the accumulation risk already present in this profile. Conclusion While no generative algorithm currently offers a perfect solution, AVATAR and Gaussian Copula represent the most viable candidates, being the only methods to satisfy both macroscopic indistinguishability and attack-based privacy criteria. These findings highlight the necessity of a structured, two-stage validation framework and suggest that, when coupled with therapeutic drug monitoring, synthetic datasets could significantly enhance multicentre collaboration while maintaining strict regulatory compliance

BackgroundThe discontinuation of Fasiglifam (TAK-875), a GPR40/FFAR1 full agonist, during Phase 3 clinical trials due to hepatotoxicity led to widespread abandonment of GPR40 as a viable therapeutic target for type 2 diabetes mellitus (T2DM). However, mechanistic evidence suggests that Fasiglifams hepatotoxicity arises from mitochondrial liability driven by high lipophilicity (aLogP = 5.31), rather than from on-target GPR40 signaling. We hypothesized that target-level failure was incorrectly inferred from compound-level safety concerns, and that superior candidates exist within publicly available databases. MethodsWe queried ChEMBL Release 36 (28 GB SQLite, 74 tables) for all compounds with documented GPR40/FFAR1 activity (UniProt: O14842). Compounds were filtered by EC50 [&le;] 10 nM in nM units with standard relation "=". Drug-likeness was assessed using Lipinskis Rule of Five (Ro5), aLogP, molecular weight (MW), hydrogen bond donors/acceptors (HBD/HBA), and polar surface area (PSA). A parallel analysis of Therapeutic Target Database (TTD v10.1.01, 4,298 targets) provided clinical context. A real-world evidence (RWE) patient stratification framework was constructed using EMR data from tens of millions of patients with >10 years of longitudinal follow-up. ResultsOf 2,637 GPR40-active compounds in ChEMBL 36, 526 (19.9%) demonstrated EC50 < 100 nM and 102 (3.9%) demonstrated EC50 < 10 nM. Eight compounds met stringent drug-likeness criteria (Ro5 violations = 0, aLogP < 5.0, EC50 [&le;] 1 nM). The lead compound (CHEMBL4859651) exhibited EC50 = 0.04 nM (8.75-fold more potent than Fasiglifam), MW = 297 Da (43% lower), and aLogP = 4.30 (19% lower), with zero Ro5 violations. Mean MW of the eight candidates was 317 {+/-} 28 Da versus 524 Da for Fasiglifam. A parallel GCK analysis identified a protein-protein interaction target (CHEMBL3885579, GCK-GKRP interface) harboring 40 exclusive compounds as an orthogonal strategy for partial GCK activation. ConclusionsSystematic cheminformatic analysis reveals that compounds with substantially superior activity and drug-likeness profiles relative to Fasiglifam exist within ChEMBL 36. Fasiglifams hepatotoxicity is attributable to compound-specific physicochemical properties, not GPR40-mediated toxicity. RWE patient stratification may further mitigate hepatotoxicity risk for next-generation GPR40 agonists. These findings argue for systematic reappraisal of GPR40 as a viable therapeutic target for T2DM.

Drug-induced seizures remain a major safety concern in drug development, yet human seizure liability is difficult to predict using conventional preclinical models. Here, we evaluated whether spontaneous calcium network activity in human induced pluripotent stem cell-derived CNS-3D Brain Organoids could predict clinically observed seizure risk across a pharmacokinetically anchored drug set. CNS-3D organoids contained neuronal and astrocytic populations, expressed neuroactive receptor and ion-channel gene programs that aligned with human cortical tissue, and exhibited reproducible spontaneous calcium oscillations across production batches. A retrospective drug panel of 66 small-molecule drugs was assembled from human clinical evidence, including 30 seizure-associated drugs and 36 comparator drugs without documented clinical seizure liability. Drugs were tested across concentration ranges anchored to reported clinical Cmax, and calcium time-series responses were integrated with chemical structure features using a machine-learning workflow. The final model predicted clinical seizure liability with an AUROC of 0.872, achieving 83.3% sensitivity and 88.9% specificity in drug-level cross-validation. Model scores also stratified seizure-associated drugs by clinical context and prevalence, suggesting that CNS-3D activity profiles capture clinically meaningful differences in seizure risk. Compared with published in vitro and preclinical seizure-liability models, CNS-3D organoid-based predictions showed improved balanced sensitivity and specificity. These findings support high-throughput calcium profiling in human CNS-3D organoids as a scalable, exposure-aware platform for predicting human seizure liability and contributing functional human data to neuro-safety assessment.

Population pharmacokinetic (popPK) models are commonly developed using ordinary differential equations (ODEs) to describe deterministic concentration-time profiles, with unexplained variability typically attributed to interindividual variability or residual error. When model misspecification is present, system-level deviations may be absorbed into these conventional variability terms, making the source and magnitude of model inadequacy difficult to assess quantitatively. Stochastic differential equations (SDEs) provide an alternative framework by introducing an explicit system-noise component into the structural model, allowing model-data mismatch to be evaluated more directly. However, historical implementation of SDE-based models in NONMEM has been technically challenging. The availability of the Fortran plug-in subroutine SDE.f90 substantially lowers this barrier and enables more practical implementation of SDE-based models in NONMEM. In this work, SDE-based nonlinear mixed-effects models were evaluated as a quantitative diagnostic framework for probing popPK model misspecification. The SDE.f90 implementation was first verified using simulated one-compartment intravenous bolus datasets with stochastic process noise. Additional simulation-estimation scenarios were then conducted under intentionally misspecified structural or stochastic assumptions, including time-varying elimination, compartmental misspecification, and residual error misspecification. Across these scenarios, the estimated system-noise parameter was generally sensitive to misspecification, with larger values usually associated with greater structural or stochastic mismatch. SDE-based modeling also helped partially separate system-level variability from residual variability and, in selected settings, supported localization of misspecification to specific model components, thereby helping guide model refinement. Overall, SDE-based popPK modeling is a useful addition to the pharmacometric diagnostic toolbox, with system-noise estimates best interpreted together with structural model evaluation, residual diagnostics, parameter behavior, and pharmacologic plausibility.

BackgroundAlkaptonuria (AKU) is a rare autosomal recessive metabolic disorder caused by deficiency of homogentisate 1,2-dioxygenase (HGD), resulting in systemic accumulation of homogentisic acid (HGA), ochronosis, and progressive multisystem disease. Although nitisinone (NTBC) lowers HGA levels, it does not correct the underlying genetic defect and induces hypertyrosinemia, highlighting the need for curative treatment approaches. We evaluated liver-directed adeno-associated virus (AAV)-mediated HGD gene therapy as a potential treatment for AKU. MethodsHgd-deficient (Hgd-/-) mice received liver-directed AAV2/8 vectors expressing codon-optimized human HGD under a liver-specific promoter. Reporter vectors were first used to assess hepatic biodistribution and transduction efficiency. Therapeutic efficacy was subsequently evaluated following AAV2/8-HGD administration (1 x 1012 vg/mouse). HGD expression was assessed by DNAscope, Western blotting, and RT-qPCR. Metabolic correction was determined using targeted LC-MS/MS and untargeted LC-HRMS metabolomics and compared with NTBC-treated Hgd-/- mice. ResultsReporter studies demonstrated liver-predominant transduction, with dose-dependent hepatocyte transduction reaching 89-93% at the highest dose. AAV2/8-HGD treatment produced robust hepatic HGD expression, with codon-optimized human HGD transcript levels approximately 33-fold higher than endogenous murine Hgd expression. Twelve weeks after treatment, plasma and urinary HGA levels were significantly reduced, with plasma HGA restored to near wild-type concentrations. Untargeted metabolomics further demonstrated marked reductions in HGA-derived phase I and II metabolites and revealed significant modulation of tricarboxylic acid cycle metabolism, consistent with partial restoration of metabolic homeostasis. Compared with NTBC-treated mice, AAV2/8-HGD achieved comparable plasma HGA reduction without elevation of upstream tyrosine pathway metabolites. ConclusionsLiver-directed AAV2/8-HGD gene therapy achieved substantial biochemical correction in Hgd-/- mice and restored metabolic flux without inducing hypertyrosinemia. These findings provide proof-of-concept supporting AAV-mediated HGD replacement as a promising long-term therapeutic strategy for AKU.

Prader-Willi syndrome (PWS) is a rare multisystemic disorder characterized by obesity, endocrine dysfunctions, and psychiatric comorbidities, which imply frequent use of psychotropic medications. They account for atypical responses to standard dosages of psychiatric drugs. Pharmacogenetics could be part of the reason for this situation, potentially offering a valuable tool for individualized treatment. This study analyzed allelic and phenotypic frequency distributions of five of the main cytochrome P450 enzymes (CYP2D6, CYP2B6, CYP2C19, CYP2C9, CYP3A4) involved in psychiatric drug metabolism in 47 patients with genetically confirmed diagnosis of PWS and compared them to reference frequencies in the general European population. Allelic frequency comparisons between the European reference population and the overall PWS cohort revealed a significant global difference for CYP2B6, with CYP2C19 and CYP2D6 showing trends toward significance. Although no global allelic differences remained significant after false discovery rate correction, post-hoc analyses consistently identified an enrichment of reduced- or non-functional alleles CYP2B619 and CYP2D610 in patients with PWS. Predicted metabolizer phenotype analyses showed a significant shift toward intermediate metabolizers of CYP3A4 in the PWS cohort, with corresponding depletion of normal metabolizers. Subgroup analyses indicated that allelic differences were more pronounced in maternal uniparental disomy and non-deletion subtypes, particularly for CYP2B6, although no significant differences were observed between PWS genetic subtypes. Overall, results imply potential differences in metabolizing activity in PWS patients, and subsequent implications in drug efficacy and tolerability. These results support the idea that pharmacogenetic testing may improve therapeutic decision-making in PWS for psychiatric treatment. Larger studies are needed to confirm these preliminary results.

Gray-market retatrutide use is increasing, but patient safety experiences remain poorly characterized. This cross-sectional analysis examined Reddit posts and comments from retatrutide-specific and broader peptide or weight-management communities through December 2025. A validated large language model classified self-reported retatrutide use and extracted author-attributed symptoms mapped to MedDRA Preferred Terms. Among 13,589 users reporting current use, 7,823 had at least one mapped symptom after exclusions. Unlike phase 2 trial findings dominated by gastrointestinal events, Reddit reports most often described appetite increase, fatigue, increased energy, nausea, food craving, insomnia, and elevated heart rate. Findings are hypothesis-generating and warrant pharmacovigilance attention.

Rationale: Efficacious dose selection for anti-tuberculosis drugs has traditionally relied on achieving plasma exposures above the minimum inhibitory concentration, but this approach has not consistently aligned with clinical outcomes. Objectives: We sought to identify early pharmacokinetic-pharmacodynamic targets most predictive of clinical efficacious dose. Methods: We conducted a back-translational, pharmacokinetic-pharmacodynamic simulation-based analysis of 15 anti-tuberculosis drugs. Using pharmacokinetic data from multiple biological matrices and a range of pharmacodynamic metrics, we established candidate exposure-response targets for attainment. We systematically evaluated the predictive accuracy of each target pair against established clinical doses to formulate a decision-making framework linking key drug properties to the most predictive targets. Measurements and Main Results: Depending on the target used, projected clinical doses varied widely - both within and across compounds - highlighting the importance of target selection for dose projection and go/no-go decisions. In general, targeting cellular lesion-level drug exposures relative to in vivo preclinical potency provided an effective approach for early dose selection. However, for highly penetrating drugs, targeting site-of-action therapeutic exposures in the caseum was more predictive of clinical dose. Based on these findings, we developed a preliminary dose prediction tool that enables drug developers to estimate clinically relevant dose ranges of compounds using in vitro and early in vivo data. Conclusions: This work establishes and validates a simple, evidence-based framework to standardize early translational decision-making on dose selection of anti-tuberculosis candidates in development.

Background: Dihydropyridine calcium channel blockers (DHP-CCB) are widely prescribed antihypertensives whose adverse effects may trigger unnecessary prescribing of additional medications, termed prescribing cascades (PC). We aimed to identify potential DHP-CCB-induced PCs using high-throughput sequence symmetry analysis (HTSSA). Methods: Using Medicare claims data (2011-2020), we identified new users aged [&ge;]66 years with continuous enrollment [&ge;]360 days before and [&ge;]180 days after DHP-CCB initiation. We screened for initiation of 446 "marker" drug classes within {+/-}90 days of DHP-CCB initiation. Sequence ratios compared marker drug initiation after versus before DHP-CCB initiation. Adjusted sequence ratios (aSR), accounting for prescribing trends over time, were calculated with 95% CIs >1 considered statistically significant. Clinical experts classified statistically significant signals as potential PCs through consensus. Results: Among 388,862 DHP-CCB initiators (mean age 76.6 {+/-} 7.5 years; 62.5% women, 92.3% with hypertension), 82 of 446 marker drug classes had significantly elevated aSRs, of which 24 were classified as potential PCs. Strongest signals ranked by highest aSR included other systemic hemostatics (aSR 2.99; 95% CI, 1.10-8.16), other nasal preparations (aSR 1.99; 95% CI, 1.47-2.70), and drugs used in erectile dysfunction (aSR 1.85; 95% CI, 1.27-2.70). Other clinically relevant signals, ranked by number needed to harm (lowest to highest), included sulfonamides (NNTH 104; 95% CI, 98-111), electrolyte solutions (NNTH 216; 95% CI, 196-241), and osmotically acting laxatives (NNTH 710; 95% CI, 540-1056). Conclusion: Potential PCs identified in this Medicare cohort reflected known and underrecognized adverse effects of DHP-CCBs. Further studies are needed to evaluate the clinical consequences of these PCs.

Background: Electrolyte replacement is ubiquitous in the acute care setting, but its familiarity cannot belie that even small dosing errors with potassium can cause lethal cardiac arrhythmias. Recently, MedAgentBench offered a benchmark for agentic artificial intelligence (AI) including the ability to correctly dose potassium based on a single rule; however, this does not adequately reflect the clinical complexity or safety concerns of an agent that has been used as the lethal injection. The purpose of this analysis was to a probe leaderboard large language model (LLM) capabilities to follow basic dosing rules to safely replace potassium in a series of clinician-annotated cases. Methods: Using a clinician panel, we developed a series of dosing principles and 20 clinical cases reflective of the complexity of potassium replacement. External clinicians were surveyed to assess practice variability and agreement to clinician panel answers. We tested GPT-5-chat with each case in triplicate, with and without the clinician curated dosing principles, and prompted the model to answer six questions involving potassium goals, dosing, route, lab frequency, concurrent interventions, and the model's perceived level of confidence for the output and complexity of the case. The primary outcome was the rate of appropriate recommendations in comparison to clinician answers. Results: A total of 54 clinicians reviewed the 20 hypokalemia cases and hypokalemia dosing guideline. Clinicians expressed "highly agree" or "somewhat agree" for 66.8% of the cases evaluated when asked if they agree with the guideline-recommended management. When given the potassium dosing guideline, total errors dropped from 165 to 104, and average accuracy improved from 45% to 65% with GPT-5-Chat. GPT-5-Chat conveyed a high level of confidence for 100% of responses, while labeling 80% and 76% of cases as highly complex with and without the criteria, respectively. Potential harm scores were considerable in both groups, however, a notable reduction in severity scores occurred with the dosing guidance document. Recommendations on concurrent interventions and dosing had the highest rate of errors in both groups. Conclusions: Benchmarks must appropriately reflect clinical complexity to be considered valuable for the deployment of agentic artificial intelligence tools in the healthcare domain. GPT-5-Chat assessment on a comprehensive medication management task for potassium replacement showed improvement with dosing guidance, yet unfit benchmarking performance.

Background: GLP-1 receptor agonists (GLP1-RAs) are an established treatment for type 2 diabetes mellitus (T2DM) and obesity. Their widespread use is set to increase through both indication expansion and patent expiry. As well as efficacy, it is crucial to understand the safety of this drug class to enable optimal use. Here we demonstrate how a genetic approach can augment signal-detection and post-market authorization surveillance. Methods: We used single nucleotide polymorphisms (SNPs) in GLP1R to recapitulate the effect of agonism with GLP1RAs on circulating glucose, glycated hemoglobin (HbA1c), body mass index (BMI) and risk of type 2 diabetes (T2DM) using Mendelian randomisation. We then tested if the adverse effect highlighted by medicines regulators of pancreatitis and the emerging effect of sarcopenia were causally related to GLP1R agonism, using this approach. Analyses were conducted in UK biobank and replicated in FinnGen and All of Us, results being combined using meta-analysis. Analyses were further stratified by a priori risk factors of age and alcohol consumption. Results: Genetically proxied GLP-1R agonism was associated with a reduction in glucose (exp({beta}) = 0.95 95% CI [0.94, 0.97]), HbA1c (exp({beta}) = 0.94 95% CI [0.92, 0.95]), and BMI (exp({beta})=0.98 95% CI [0.97, 0.99]); and a reduced risk of T2DM (OR = 0.82 95% CI [0.79 to 0.86]). Risk of acute and chronic pancreatitis was however increased (OR = 1.10 95% CI [1.01 to 1.20] and OR = 1.05 95% CI [0.95, 1.17], respectively), which varied as a function of age with risk most pronounced in those aged 50-59 years-old (OR = 1.79 95% CI [1.43, 2.24], OR = 1.57 95% CI [1.16, 2.12]) and in drinkers (OR = 1.32 95% CI [1.12, 1.54], OR = 1.36 95% CI [1.12, 1.65]). Risk of sarcopenia also increased (OR 1.34; 95% CI 1.05,1,71). Conclusions: Genetically proxied agonism with GLP-1RAs recapitulated the pharmacological effects of GLP1-1RAs on glycaemic traits, BMI and T2DM risk. This approach supports a causal effect of GLP-1RAs on the well reported adverse effects of pancreatitis and further indicates age and alcohol consumption as risk modifying effects. The less well reported but emerging effect of sarcopenia appears to also be casually related to agonism at GLP-1R. These analyses suggest a genetic approach could be used as an adjunct to signal detection studies to enhance safety regulation as well as personalisation of the use of these drugs.

Clinical pharmacokinetic (PK) modelling is constrained by sparse sampling, limited general-isability of single-drug models, and labour-intensive workflows, making it difficult to infer complete drug exposure from limited concentration observations. We present the Pharmacokinetic Foundation Model (PKFM), a grey-box Transformer framework pre-trained across 32 drugs that reconstructs concentration-time profiles from sparse concentration observations, dosing events, molecular descriptors, and physiological covariates while preserving output interpretability. In representative oral PK curves, three sparse input points recovered the principal absorption-elimination trajectory, achieving coefficient of determination (R2) = 0.992 for Midazolam oral and R2 = 0.990 for Verapamil oral. Using reconstructed curves in NONMEM (nonlinear mixed-effects modelling) improved covariance stability and individual prediction accuracy. Contrastive-learning embeddings supported Top-10 physiologically based pharmacokinetic (PBPK) candidate retrieval, with 75.6% of observations within the 2-fold range. A pharmacometrics-informed AI Agent (PM Agent) outperformed general-purpose programming tools in stability and pairwise win rate on a standardised modelling benchmark, with each run requiring human pharmaco-metrician confirmation before downstream use. These results support cross-drug pre-trained PK models as an information-completion layer for sparse PK evidence and a structured scaffold for the modelling workflow; clinical or regulatory use requires prospective validation, broader external benchmarking, and independent expert assessment.

Introduction: Acute low back pain is a leading cause of disability worldwide. Clinical guidelines recommend non-pharmacological therapies as first-line treatment and advise caution with opioid prescribing. However pharmacological therapies, including opioids and gabapentinoids, remain commonly used. The comparative risks of subsequent opioid use disorder (OUD) and overdose diagnosis associated with initial treatment modality in large, real-world populations is not well characterized. We estimated the incidence of new-onset OUD and overdose diagnosis among opioid-naive, Medicaid-insured adults with newly diagnosed acute low back pain and estimated the association between initial treatment modalities and subsequent OUD and overdose diagnosis risk. Methods: We conducted a retrospective cohort study using Medicaid T-MSIS Analytic files from 25 states (2016-2019). We identified opioid-naive adults with a new diagnosis of acute low back pain who initiated pharmacologic or non-pharmacologic treatment within 1 month of diagnosis. The primary outcome was incident OUD and overdose diagnosis (based on diagnosis codes in claims) during follow-up. Associations between initial treatment modality and OUD and overdose diagnosis risk were estimated using a non-parametric, doubly robust estimator to adjust for measured confounding. Results: The cohort included 525,002 opioid-naive adults initiating treatment for low back pain. The cumulative incidence of OUD and overdose diagnosis was 1.5% and 2.4% at 7 and 13 months, respectively. Compared to non-use, use of gabapentinoids during the first month of treatment was associated with the highest relative risk (increasing risk) by 130.1%, 95% confidence interval (CI): 117.8%, 142.3%), the second-highest relative risk was estimated for higher-dose opioids, defined as > 50 daily Morphine Milligram Equivalents (MME) (118.1%, 95% CI: 99.2%, 137.0%). Lower-dose, short-duration opioids ([&le;] 50 MME, [&le;] 7 days) were also associated with elevated risk, though substantially smaller in magnitude (20.8%, 95% CI: 13.8%, 27.9%). In contrast, non-pharmacologic, non-interventional therapies were associated with reduced OUD and overdose diagnosis risk, with physical therapy demonstrating the largest relative reduction of 34.0% (95% CI: -40.9%, -27.1%). Discussion: In opioid-naive Medicaid patients with acute low back pain, initial non-pharmacologic treatment was associated with reduced OUD and overdose diagnosis risk. Gabapentinoids and opioids were each associated with increased risk; for opioids, the degree of risk increased with higher doses and durations. These results support guideline recommendations favoring non-pharmacologic treatment as first-line therapy and indicate the importance of cautious prescribing when pharmacologic treatment is considered.

Background: Psychedelics are emerging as potential management options for chronic musculoskeletal pain due to preliminary evidence of effectiveness and low addictive potential, but patients perceptions remain unknown. This study assessed patient perceptions regarding psilocybin for musculoskeletal pain. Methods: We conducted a cross-sectional survey of adults ([&ge;]19) with musculoskeletal pain attending a hospital-based orthopaedic clinic. Participants reported demographics, perceptions of psychedelics for pain management, and willingness to participate in psychedelic research. Multivariable regression explored factors associated with perceived analgesic potential, and willingness to try a full therapeutic dose of psilocybin or a microdose. Results: Among 295 participants, 73% reported moderate-to-severe pain; 75% used analgesics; of these, 41% used opioids (86/209). While 24% reported prior psychedelic use, only 3% had discussed psychedelics with a healthcare provider. Most perceived that psilocybin had moderate-to-high effectiveness for pain (76%). Most respondents endorsed a moderate-to-high willingness to try microdoses (58%) and macrodoses (53%) of psilocybin for pain. Prior non-therapeutic psychedelic use predicted a 1.05-unit increase in perceived analgesic potential on the 10-point scale (p=.013). Willingness to try a macrodose of psilocybin was most strongly associated with prior non-therapeutic (B=3.16) and therapeutic (B=2.42) psychedelic use; in contrast, pain severity had a significant but modest association, with a 0.21-point increase in willingness for every 1-unit increase in pain severity (p=.017). Similarly, willingness to try a microdose of psilocybin was predicted by non-therapeutic (B=2.82) and therapeutic (B=2.48) use, whereas the effects of pain severity (B=0.20) and younger age (B=-0.30) were significant but small. Most respondents (52%) reported moderate-to-high willingness to participate in a trial of psilocybin for pain relief, and health risks were the primary concern (33%). Conclusions: Study findings suggest a majority hold neutral-to-positive perceptions of psilocybin for pain. Addressing perceived barriers, including health effects and gaps in patient knowledge, should be considered when designing future trials.

Background and Aims Hypertension during pregnancy is a major cause of maternal and neonatal morbidity and mortality, yet the efficacy and safety of antihypertensive treatments in this setting remain uncertain. We evaluated the effects of antihypertensive drug targets on adverse pregnancy-related outcomes using genetic variants to instrument target perturbation. Methods We performed drug target Mendelian randomization to mimic pharmacological perturbation of targets from six commonly used antihypertensive drug classes, using data from up to 671,922 pregnant women. Genetic variants near drug target genes associated with systolic or diastolic blood pressure were selected as instruments. We estimated effects of target modulation on six primary and eight secondary pregnancy outcomes. Results Genetically instrumented downregulation of blood pressure through beta-blocker (BB) and calcium-channel blocker (CCB) targets, particularly ADRB1 and CACNB2, was associated with a reduced risk of hypertensive disorders of pregnancy, including preeclampsia. For example, CACNB2-instrumented lowering corresponded to a 7% (95% CI: 5-9%) reduction in preeclampsia risk per 1 mmHg decrease in blood pressure. For most other targets, estimates were directionally consistent but imprecise. Across additional outcomes, effects varied by target, with suggestive evidence for reduced risks of miscarriage, preterm birth, small-for-gestational-age birth, and labour induction, although these estimates were accompanied by substantial uncertainty. Conclusions These findings support a protective effect of BB and CCB targets on hypertensive disorders of pregnancy and highlight potential target-specific differences in safety. This work illustrates the value of Mendelian randomization in addressing clinical uncertainties where robust trial evidence is limited.

NAT2 is an important pharmacogene which encodes the N-acetyltransferase 2 enzyme that is involved in the metabolism of multiple medications, and variants in this gene can affect patient response to these medications. CPIC has published a clinical guideline for prescribing hydralazine using NAT2 genotypes. Just prior to the guideline, updated NAT2 star allele numbering and definitions were released, differing somewhat from the historical nomenclature. Clinical pharmacogenomic testing panels often test for the most common star alleles, so knowledge of the most common updated NAT2 star alleles is critical for the implementation of the CPIC NAT2/hydralazine guideline. We first determine NAT2 diplotype frequencies from UK Biobank (UKBB) 200k phased genomes, then analyzed allele, diplotype, and phenotype population frequencies from the All of Us Research program, PennMedicine BioBank (PMBB) and UKBB 500k datasets. We found that analyzing NAT2 diplotypes from phased data provides critical information for algorithms designed to predict diplotypes from unphased data. We observed that NAT2*5, *6, and *4 were the most common star alleles in that order, and the top 11 most frequent NAT2 star alleles were the same across all biobanks. However, differences in star allele frequencies across biogeographical populations were observed. The largest difference led to a higher frequency of NAT2 poor metabolizer phenotypes as compared to rapid and intermediate metabolizer phenotypes in all global populations except in the EAS population, where NAT2 poor metabolizers were in the minority.

BackgroundPrioritizing therapeutics from transcriptomic data remains a key challenge in precision medicine. Signature reversal approaches, most commonly implemented through Gene Set Enrichment Analysis (GSEA), have been widely used to match disease signatures to candidate drugs. However, enrichment-based methods can be sensitive to noise and are restricted to previously profiled compounds MethodsWe developed RANKOR, a machine-learning framework designed to rank candidate drugs directly from transcriptomic signatures. Rather than predicting full expression profiles, RANKOR learns structured latent representations of transcriptional responses alongside chemical structure, enabling prioritization from standardized signatures derived from disease states or treatment perturbations. The framework is applicable to both bulk and single-cell transcriptomic data. ResultsAcross large-scale perturbational datasets, RANKOR achieved consistently lower median ranks than similarity- and distance-based approaches, while showing performance comparable to, and in some settings improved over, GSEA. The model generalized across unseen cell types and retained performance in single-cell settings, where it provided more consistent prioritization than existing approaches, such as ASGARD. RANKOR further enabled prioritization of transcriptionally unseen compounds through chemical-space embedding and achieved substantially reduced computation times. Robustness analyses demonstrated stable performance under moderate noise and degradation under extreme perturbation or gene shuffling. Gene attribution analyses indicated that prioritization decisions are driven by coherent and mechanism-relevant transcriptional programs. ConclusionsRANKOR provides a scalable framework for transcriptomics-guided drug prioritization that can complement and extend existing approaches, such as GSEA. It can also support therapeutic hypothesis generation from bulk and single-cell data while leveraging the generalisability and computational efficiency of machine learning models.